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Abstract 

A study of the integrability of one-dimensional quantum mechanical 
many-body systems with general point interactions and boundary con- 
ditions describing the interactions which can be independent or depen- 
dent on the spin states of the particles is presented. The corresponding 
Bethe ansatz solutions, bound states and scattering matrices are explic- 
itly given. Hamilton operators corresponding to special spin dependent 
boundary conditions are discussed. 

Exactly solvable models of a single quantum particle moving in a local singular poten- 
tial concentrated at one or a discrete number of points have been extensively discussed in 
the literature, see e.g. [[I], |2|, [| and references therein. In the one dimensional case, the 
local singular potential (contact interactions) at, say, the origin (x = 0) can be character- 
ized by the boundary conditions imposed on the wave function (p at x = 0. There are two 
classes of such boundary conditions: separated and nonseparated boundary conditions, 
corresponding to the cases where the perturbed operator is equal to the orthogonal sum 
of two self-adjoint operators in L 2 (— oo,0] and L 2 [0,oo) and when this representation is 
impossible, respectively. The many-body problems with pairwise interactions given by 
such boundary conditions are generally not exactly solvable [|], || |J. In the present pa- 
per we give a systematic description for integrable models of many-body systems with 
pairwise interactions given by such singular potentials for the case where the boundary 



conditions are independent as well as for the case where they are dependent on the spin 
states of the particles. 

We first consider the case of spin independent boundary conditions. The family of 
point interactions for the one dimensional Schrodinger operator — can be described by 
unitary 2x2 matrices via von Neumann formulas for self-adjoint extensions of symmetric 
operators, since the second derivative operator restricted to the domain Qj°(R\ {0}) has 
deficiency indices (2, 2). The nonseparated boundary conditions describing the self-adjoint 
extensions have the following form 

where 

ad — be = 1, 9,a,b,c,d G M. (2) 

<p(x) is the scalar wave function of two spinless particles with relative coordinate x. (jl]) 
also describes two particles with spin s but without any spin coupling between the particles 
when they meet (i.e. for x = 0), in this case ip represents any one of the components 
of the wave function. The values 6 = b = 0, a = d = 1 in ([!]) correspond to the case 
of a positive (resp. negative) 5-function potential for c > (resp. c < 0). For general 
a, b, c and d, the properties of the corresponding Hamiltonian systems have been studied 
in detail, see e.g. |7], [J. 

The separated boundary conditions are described by 

<f/(0+) = q + <p(0+) , p'(0_) = g>(0_), (3) 

where q ± G M U {oo}. q + = oo or q~ = oo correspond to Dirichlet boundary conditions 
and q + = or q~ = correspond to Neumann boundary conditions. 

To study the integrability of one dimensional systems of iV-identical particles with 
general contact interactions described by the boundary conditions (|I|) or (|3|) that are 
imposed on the relative coordinates of the particles, we first consider the case of two 
particles (N = 2) with coordinates x±, Xi and momenta ki, k 2 respectively. Each particle 
has n-'spin' states designated by s± and S2, 1 < Sj < n. For x± ^ x 2 , these two particles 
are free. The wave functions tp are symmetric (resp. antisymmetric) with respect to the 
interchange (xi,Si) <-> (2:2, Sa) f° r bosons (resp. fermions). In the region X\ < x 2 , from 
the Bethe ansatz the wave function is of the form, 

ip = Ul2e i(kix 1+ k 2 x 2 ) + U2ie i(k 2 x 1+kl x 2 ) ^ 

where u\ 2 and u 2 \ are n 2 x 1 column matrices. In the region x% > x 2 , the wave function 
has the form 

<p = (P 12 Ul2 )e i{klX2+k2Xl) + (P 12 u 21 )e i{k2X2+klXl \ (5) 
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where according to the symmetry or antisymmetry conditions, P 12 = p 12 for bosons and 
P 12 = —p 12 for fermions, p 12 being the operator on the n 2 x 1 column that interchanges 
Si <-> s 2 . 

Let ki2 = (fcx — fc 2 )/2. In the center of mass coordinate X = {x\ + x 2 )/2 and the 
relative coordinate x — £ 2 — x%, we get, by substituting and (|) into the boundary 
conditions at x — 0, 

«21 = Y^Un , (6) 

12 = 2ie ie k 12 P 12 + »fc 12 (q - d) + (A:i 2 ) 2 & + c 

21 ^ 12 (a + d) + (A; 12 ) 2 6 - c U 



for boundary condition and 



v i2 _ ifeia + 9 
r 21 



for boundary condition @, where q = q + = —q_ £ MU {00}. 

For N > 3 and Xi < £2 < ... < x^, the wave function is given by 

(9) 

+ (N\ — 2) other terms. 

The columns u have x 1 dimensions. The wave functions in the other regions are 
determined from (|S|) by the requirement of symmetry (for bosons) or antisymmetry (for 
fermions). Along any plane Xi = Xj+i, i 6 1,2, ...,N — 1, from similar considerations as 
above we have 

OLa-LOLi-.-onoti+i— ajv — ^a i+iai M ai02---Oi+iOi---ajv ) (10) 

where 

yii+i _ 2ie k ai g i+1 P + ik aiai+1 {a — d) + (fca jQ ^ +1 ) 6 + c 
ik aiai+1 (a + d) + (k aiCti+1 ) 2 b — c 
for nonseparated boundary condition and 

yii+l _ + , > 



for separated boundary condition. Here = (k ai — k cti+1 )/2 play the role of spectral 

parameters. p n+1 = f or bosons and P u+1 = —p u + l f or fermions, with p u+1 the 

operator on the n N x 1 column that interchanges Sj <-> Sj+i. 

For consistency K must satisfy the Yang-Baxter equation with spectral parameter 
ig, 0,0,^,0, i.e., 



-i/-m,m+l-\rm+l,m+2-\rm,m+l -i^m+l,m+2-i^m,m+l-i^m+l,m+2 

ij 1 kj ki ki kj ij > 

or 

r ij r fcj — x ki 1 kj I ij 
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if m, r, s are all unequal, resp. 



^3 3 ^ ' ^3 ^.7 

if m, r, s, g are all unequal. 

The operators Y given by (|TTD satisfy the relation flHD for all a, b, c, d However the 
relations ( |13D are satisfied only when 9 = 0, a = d and 6 = 0, that is, according to the 
constraint (Q), 9 — 0, a = d — ±1, 6 = 0, c arbitrary. The case a = d = 1, 9 = b = 



corresponds to the usual ^-function interactions, which has been investigated in [1C, 11 



T2|j . The case a — d — — 1, 9 = b = is related to another singular interactions between 
any pair of particles (for a — d — — 1 and 9 = b = c = see (7|, ||), which is in fact unitarily 
equivalent to the ^-interaction, under a non-smooth "kink type" gauge transformation 
U = I\i>j sgn(a;j — Xj). Associated with the separated boundary condition, the operators 
Y given by (|l~2|) satisfy both the relations ( |i~3l) and ([14]) for arbitrary q. Therefore with 
respect to iV-particle (either boson or fermion) problems, there are two non-equivelant 
integrable one parameter families with contact interactions described respectively by one 
of the following conditions on the wave function along the plane Xi = Xj for any pair of 
particles with coordinates and Xj, 

y?(0+) = +p(0_), <p'(0+) = c^(0_) + p'(0_) , c e R ; (15) 

<f/(0+) = q<p(0+), p'(0_) = -q<p(0-) , q E RU {oo} . (16) 

The wave functions are given by (§) with the it's determined by (|T0| ) and initial conditions. 
The operators K in (|H]) are given respectively by 



yii+l __ K^on kg t+1 )P ll+1 + C _ 

and 



yii+l _ ^Qj+l) + ^ 

1+1 1 H^ai — k ai+1 ) — 2q 
When q < 0, there exist 2 N<yN ~ 1 ^ 2 bound states for the case (18) of separated boundary 



conditions, with wavefunction 

^N,e = u e _l[(9(x k - xi) + e kl 9(x t - x k ))e q ^>i lx '~ x ^ (19) 

k>l 

and eigenvalue E = —q 2 N(N 2 — l)/3, where u e is the spin wave function and e = 
{e kl : k>l}; e kl = ±, labels the 2 N ^ N ~ 1 ^ 2 -Md de generacy. 

We consider now the case of spin dependent boundary conditions. For a particle with 
spin s, the wave function has n — 2s + 1 components. Therefore two particles with contact 
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interactions have a general boundary condition described in the center of mass coordinate 
system by: 

f if) \ ( A B \ ( i\) \ 



v )o + \ c D J V ^ lo- 



rn 



where ip and ip' are n 2 -dimensional column vectors, A, B, C and D are n 2 x n 2 matrices. 
The boundary condition ( PUD can include not only the usual contact interaction between 
the particles, but also a spin coupling of the two particles if the matrices A, B, C, D are 
not diagonal. 

The matrices A, B, C, and D are subject to restrictions due to the required symmetry 

d 2 

condition of the Schrodinger operator. For any u, v G C°°(-K\{0}), < —-j-^u, v >L 2 (M,<u n ) 
d 2 

— < u, — — v >l 2 (]R (V n )= 0, which, together with (pO) imply 

dx 2 ' 

A ] D-C ] B=l, B^D = D ] B, A j C = C^A, (21) 

where f stands for the conjugate and transpose. Obviously (|1|) is the special case of ( |20| ) 
for s = 0. 

In the following we study quantum systems with contact interactions described by 
the boundary condition (]2T)|), in particular, iV-body systems with S- interactions. We first 
consider two spin-s particles with ^-interactions. The Hamiltonian is then of the form 

^ =( "S~S )l2 + 2M(Xl_X2) ' (22) 

where I2 is the n 2 x n 2 identity matrix, h is an n 2 x n 2 Hermitian matrix. If the matrix 
h is proportional to the unit matrix I2, then H is reduced to the usual two-particle 
Hamiltonian with contact interactions but no spin coupling. 

Let e a , a = 1, ...,n, be the basis (column) vector with the a-th component as 1 and 
the rest components 0. The wave function of the system (|2"2"|) is of the form 



7 P= <j>ap{x\,x 2 )e ol ®e !3 . (23) 

a,P=l 



In the center of mass coordinate system, the operator (|22|) has the form 

The functions <ft = 4>{x, X) from the domain of this operator satisfy the following boundary 
condition at x = 0, 

n 

0^(0+ X)- 0^(0- X)= £ ^(0,4 0^(O + ,X)=^(O-,X), (25) 

a,P=l 
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01,(3 = 1, n, where the indices of the matrix h are arranged as 11, 12, In; 21, 22, 2n; 
nl,n2, ...,nn. (^5|) is a special case of (|20|) for A = D = I 2 , I? = and C = h. h acts 

n 

on the basis vector of particles 1 and 2 by he a ® ep — 2_j hap^xe-y ® e\. 

7,A=1 

The wave functions are still of the forms (f|) (resp. |5|) in the region x\ < x 2 (resp. 
X\ > x 2 ). Substituting them into the boundary conditions (^5|), we get 

f U 12 + U 2 1 = P 12 (U 12 + U 21 ), 

I ikn(u 2 i - u 12 ) = hP 12 (u 12 + u 21 ) + ih 2 P 12 (u 12 - u 21 ). 



(26) 



Eliminating the term P l2 ui 2 from (|26"D we obtain the same relation as (0), u 2 i = Y^ 2 ^- 
Nevertheless the Y operator is given by 



Y 2 12 = [2ik l2 - h]- l [2ik l2 P 12 + h\. 



(27) 



For a system of N identical particles with ^-interactions, the Hamiltonian is given by 



N 



N Q2 
i=l ux i 



(28) 



where Ijv is the n x n identity matrix, hij is an operator acting on the i-th and j-th bases 
as h and the rest as identity, e.g., hi 2 = h <g> 1 3 <g> ...ljv, with lj the n x n identity matrix 
acting on the i-th basis. The wave function in a given region, say x\ < x 2 < ... < x^, is 
of the form (Pf), with 



and 



v = Y jj+1 v 

n, [^ik a]aj+1 — hjj + i] 1 [2ik Qjaj+1 P j:i+1 + hjj + i] 



+1 



ia 2 ...otj +1 aj...a N 



(29) 



(30) 



From the Yang-Baxter equations it is straightforward to show that the operator Y 
given by (^) satisfies all the Yang-Baxter relations if 



[ha, Fi- 



ts . 



(31) 



Therefore if the Hamiltonian operators for the spin coupling commute with the spin 
permutation operator, the iV-body quantum system (^3) can be exactly solved. The 

N 

wave function is then given by @ and (|29|) with the energy E — kf . 



i=l 



For the case of spin-|, a Hermitian matrix satisfying fl3~T|) is generally of the form 



/a t\ t\ c \ 

e*i f 9 e 2 

e* 9 f e 2 

V c* e* e* b J 



(32) 
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where a, b, c, /,ei,e 2 6 €, 5 6 1R. We recall that for a complex vector space V, a matrix 
R taking values in End c (V ® V) is called a solution of the Yang-Baxter equation without 
spectral parameters, if it satisfies R12R13R23 = R23R13R12, where Rij denotes the matrix 
on the complex vector space V <8> V <S> V, acting as R on the z-th and the j-th components 
and as identity on the other components. When V is a two dimensional complex space, the 
solutions of the Yang-Baxter equation include the ones such as R q which gives rise to the 
quantum algebra SU q (2) and the integrable Heisenberg spin-| chain models such as the 
XXZ model (R corresponds to the spin coupling operator between the nearest neighbor 



spins in Heisenberg spin chain models) fl3|, 14, O, 16]. Nevertheless in general does 



111 

2 

12- 

But fl3"2| ) includes the Yang-Baxter solutions, such as R q , that gives integrable spin chain 
models (for an extensive investigation of the Yang-Baxter solutions see [T7], TE\). Therefore 



not satisfy the Yang-Baxter equation without spectral parameters: hl 2 hl 3 h 23 7^ h 23 hl 3 h 2 



for an iV-body system to be integrable, the spin coupling in the contact interaction ( 281 ) 
is allowed to be more general than the spin coupling in a Heisenberg spin chain model 
with nearest neighbors interactions. 

For N = 2, from ( Effi) the bound states have the form, 

r a = u a e^^\, a = l,...,n 2 , (33) 



where u a is the common a-th eigenvector of h and P 12 , with eigenvalue A , s.t. hu a = 
A a u a and c + aA a < 0, P l2 u a = u a . The eigenvalue of the Hamitonian H corresponding 
to the bound state (|3~3"D is — (c + aA a ) 2 /2. We remark that, whereas for the case of the 
boundary condition ([[]), for a 5 interaction one has a unique bound state, here we have 
n 2 bound states. By generalization we get the bound state for the iV-particle system, 

ipa =v a e~ £± ^ SL ^ ]Xi ~ Xil , a = l,...,n 2 , (34) 

where v a is the wave function of the spin part satisfying P^v a = v a and hijV a = A a v a , 
for any % ^ j. 

It is worth mentioning that tp^ i s °f the form @ in each of the above regions. For 
instance comparing ip 1 ^ with @ in the region x\ < x 2 ... < x^ we get 

fcj = —i _ ltiL_^ (j\r _ l) ; k 2 = h + ic, k 3 = k 2 + ic, k N = -k x , (35) 

for a = 1, ...,n 2 . The energy of the bound state ip^ * s 

E a = - {c+ ^ a)2 N(N 2 -l). (36) 

Now we pass to the scattering matrix. For real k\ < k 2 < ...k^, in each coordi- 
nate region such as ii < x 2 < ...a; at, the following term in @ describes an outgo- 
ing wave ip out = U\ 2 _ N e^ klX1+ '" +kNXN \ An incoming wave with the same exponential 
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as iJj out is given by ip in = [p 1N p 2 ( N 1 \..]u N ( N ^ 1 y^ 1 e t( - kNXN+ - +klXl ^ in the region x N < 
xm-i < ■■■ < X\. From ( p9|) the scattering matrix S defined by ip out = Sip in is given by 
S = [X 21 X 31 ...X N1 ][X 32 X A2 ...X N2 }...[X N{N _ 1) }, where X l3 = Y$pv. 

The scattering matrix S is unitary and symmetric due to the time reversal invariance 
of the interactions. < s[s 2 ...s' n \S\sis 2 ...sn > stands for the S matrix element of the 
process from the state (kiSi, k 2 s 2 , k N S]y) to the state (kis[, k 2 s' 2 , k N s' N ). 

The scattering of clusters (bound states) can be discussed in a similar way as in [ p"2|| . 
For instance for the scattering of a bound state of two particles (xi < x 2 ) on a bound state 
of three particles (£3 < X4 < x$), the scattering matrix is 5* = [A32X42 A52] [X31X41X51] . 

The integrability of many particles systems with contact spin coupling interactions 
governed by separated boundary conditions can also be studied. Instead of ( PU|) we need 
to deal with the case 

0'(O+) = G+0(O + ), 0'(O_) = G-0(O_), (37) 

where are Hermitian matrices. For G + = G~ = G, G' = G, there is a Bethe Ansatz 
solution to © with Y^^ a . in (|29"D given by 



yii+l ik aiCH+1 + G , 



Let T be the set of n 2 eigenvalues of G. For any A a G T such that A a < 0, there are 
2-2V(iv-i)/2 Dounc i states for the iV-particle system, 

C = - xi) + £ki0(xi - x k ))e Xa ^»i lxi ~^, (39) 

k>l 

where v ae is the spin wave function and e = {e^ : k > Z}; = ±, labels the 2 Ar ( Ar ~ 1 )/ 2 - 
fold degeneracy. The spin wave function v here satisfies -P^tw = CijV a e for any 2 7^ j, that 
is, p l ^v ae = €ijV a e for bosons and p^tw = —€ijV ae for fermions. 

Again ip^ e is of the form (H) in each of the regions x ix < x i2 < ... < x iN . For 
instance comparing ip^ e with (||) in the region X\ < x 2 ... < x^ we get k\ = i\ a (N — 1), 
k 2 = ki — 2i\ a , k 3 = k 2 — 2i\ a ,...,k]y = —hi. The energy of the bound state ip^ e is 
E a = -X 2 a N(N 2 - l)/3. 

We have investigated the integrable models of iV-body systems with contact spin cou- 
pling interactions. Without taking into account the spin coupling, the boundary condition 
(P is characterized by four parameters (separated boundary conditions are a special lim- 
iting case of these). Obviously the general boundary condition fl20|) we considered in this 
article has much more parameters. A complete classification of the dynamic operators 
associated with different parameter regions remains to be done. As we have seen, the 



S 



case A = D = I 2 , B = 0, C = h corresponds to a Hamiltonian with 5- interactions of the 
form (|22|) (for N = 2). It can be further shown that (for N = 2) the following boundary 
condition 

f tp \ ( I B \ ( ip \ 



(40) 



corresponds to a Hamiltonian if of the form: H = —D%.(1 + BS) — BD X S', where B is an 
n 2 x n 2 Hermitian matrix, D x is defined by (D x f)(ip) = -f for / G C °°(iR/{0}) 

and (p a test function with a possible discontinuity at the origin. The boundary condition 

(41) 

V v /o- 

describes the Hamiltonian if = — -D 2 + iB(2D x 5 — 5'). 
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